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We present measurements of the B — * rf K branching fractions; for B + — > n'K + we measure 
also the time-integrated charge asymmetry A c h, and for B° — * r/'K® the time dependent CP- 
violation parameters S and C. The data sample corresponds to 232 million BB pairs produced 
by e+e" annihilation at the T(45). The results are B(B + r) K + ) = (68.9 ± 2.0 ± 3.2) x 10" 6 , 
B{B° -> -n'K ) = (67.4 ±3.3 ±3.2) x 10~ 6 , A ch = 0.033 ±0.028 ±0.005, 5" = 0.30 ±0.14 ±0.02, and 
C — —0.21 ± 0.10 ± 0.02, where the first error quoted is statistical and the second systematic. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



Measurements of time-dependent CP asymmetries 
in B° meson decays through a Cabibbo-Kobayashi- 
Maskawa (CKM) favored b — > ccs amplitude pj have 
provided a crucial test of the mechanism of CP viola- 
tion in the Standard Model (SM) 0. Such decays to a 
charmonium state plus a K° meson are dominated by 
a single weak phase. Decays of B° mesons to charm- 
less hadronic final states, such as (f>K°, K + K~ K° , ij'K , 
7T°K° and /o^SO)!^ , proceed mostly via a single pen- 
guin (loop) amplitude with the same weak phase 0] , but 
CKM-suppressed amplitudes and multiple particles in 
the loop introduce other weak phases whose contribution 
is not negligible 0, S, E H . 

For the decay B° — > r/'K , these additional contribu- 
tions are expected to be small, so the time-dependent 
asymmetry measurement for this decay provides an ap- 
proximate measurement of sin2/3. Theoretical bounds 
for the small deviation AS between the time-dependent 
CP- violating parameter measured in this decay and in 
the charmonium K° decays have been calculated with an 
SU(3) analysis 0,0- Such bounds have been improved 
by measurements of B° decays to a pair of neutral light 
pseudoscalar mesons 0, . From these and other mea- 
surements, improved model-independent bounds have 
been derived jfj , with the conclusion that AS is expected 
to be less than 0.10 (with a theoretical uncertainty less 
than ~0.03). Specific model calculations conclude that 
AS* is even smaller 7]. A significantly larger AS* could 
arise from non-SM amplitudes 0. 

The time-dependent CP-violating asymmetry in the 
decay B° — > r/' K° has been measured previously by the 
BABAR [HI and Belle £3 experiments. In this Let- 
ter we update our previous measurements with an im- 
proved analysis and a data sample four times larger . We 
also measure the B° — > r/'K and B + — » rfK + branching 
fractions 0], and for B + — ► r/ / K + the time-integrated 
charge asymmetry A c h = (r _ — r + )/(r~ ± T + ) where 
= T(B ± -> r/'K^). In the SM A ch is expected to be 
small; a non-zero value would signal direct CP violation 
in this channel. 

The data were collected with the BABAR detector Q 
at the PEP-II asymmetric e + e~ collider [l^. An in- 
tegrated luminosity of 211 fb _1 , corresponding to 232 
million BB pairs, was recorded at the T(4S) resonance 
(center-of-mass energy y/s = 10.58 GeV). Charged par- 



ticles are detected and their momenta measured by the 
combination of a silicon vertex tracker (SVT), consisting 
of five layers of double-sided detectors, and a 40-layer 
central drift chamber, both operating in the 1.5 T mag- 
netic field of a solenoid. Charged-particle identification 
(PID) is provided by the average energy loss in the track- 
ing devices and by an internally reflecting ring-imaging 
Cherenkov detector (DIRC) covering the central region. 
Photons and electrons are detected by a CsI(Tl) electro- 
magnetic calorimeter. 

From a candidate BB pair we reconstruct a B° de- 
caying into the flavor eigenstate / = r/'K° (Bcp)- We 
also reconstruct the vertex of the other B meson (B tag ) 
and identify its flavor. The difference A< = tcp — itag of 
the proper decay times tcp and £ t ag of the CP and tag 
B mesons, respectively, is obtained from the measured 
distance between the Bqp and B tag decay vertices and 
from the boost ((3j — 0.56) of the e + e~ system. The At 
distribution is given by: 

e -|At|/r 

F(At) = [l±Aw± (1) 

4t 

(1 - 2w) (S sin(Am d At) - C cos(Am d At))}. 

The upper (lower) sign denotes a decay accompanied by 
a B° (B°) tag, r is the mean B° lifetime, Arrid is the 
mixing frequency, and the mistag parameters w and Aw 
are the average and difference, respectively, of the prob- 
abilities that a true B° is incorrectly tagged as a B° or 
vice versa. The tagging algorithm 16] has seven mutually 
exclusive tagging categories of differing response purities 
(including one for untagged events that we retain for yield 
determinations). The measured analyzing power, defined 
as efficiency times (1 — 2w) 2 summed over all categories, 
is (30.5 ± 0.6)%, as determined from a large sample of 
B-decays to fully reconstructed flavor eigenstates (Bfl av ). 
The parameter C measures direct CP violation. If C = 0, 
then S = sin2/? ± AS. 

Monte Carlo (MC) simulations JjJ3| of the signal decay 
modes, BB backgrounds, and detector response are used 
to tailor the event selection criteria. We reconstruct B 
meson candidates by combining a K + or a K® with an 
r/' meson. We reconstruct rj' mesons through the decays 
i P°7 iv'ry) and i -> rrn+ir- with 77 -> 77 (7^ (rr)7r7r ) 
or i] — > 7T + 7r~7r (^(3 7r ) 7rir )- For the K + track we require 
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an associated DIRC Cherenkov angle between —5 and 
+2 standard deviations (a) from the expected value for 
a kaon. We select K® — ► tt + tt~ decays by requiring 
the 7r + 7r - invariant mass to be within 12 MeV of the 
nominal K° mass and by requiring a flight length with 
significance >3ct. We select K° — ► 7r°7r° decays by re- 
quiring that the 7r°7r° invariant mass be within 30 MeV 
of the nominal K° mass. Daughter pions from rf de- 
cays are required to have PID information inconsistent 
with proton, electron and kaon hypotheses. The photon 
energy E 1 must be greater than 30 MeV for tt° candi- 
dates, 50 (100) MeV for i] candidates for the rf t \^K Q 
( 7 7^( 77 ) 7r7r ^ + ) samples, an d greater than 100 MeV for 
r)' pl candidates. We make the following requirements 
on the invariant mass (in MeV): 490 < to 77 < 600 for 
■q — ► 77, 120 < m 77 < 150 for tt° (100 < m 77 < 155 in 
K° s -► vr 7r ), 510 < < 1000 for p°, 520 < 77W < 
570 for n vr+TT-vr , 945 < ny < 970 for ^ (rr)mr , and 
930 < m,/ < 980 for ^ r 

A _B meson candidate is characterized kine- 
matically by the energy-substituted mass ibes = 

,J (^s + p • pb) 2 /Eq — p 2 B and the energy difference 

AE = Eg — 7[y/s, where (E ,p ) and (E b ,Pb) are 
four-momenta of the T(4S) and the B candidate, respec- 
tively, and the asterisk denotes the T(45) rest frame. We 
require \AE\ < 0.2 GeV and 5.25 < m ES < 5.29 GeV. 

To reject the dominant background from continuum 
e + e~ — > qq events (q = u, d, s, c), we use the angle 0t be- 
tween the thrust axis of the B candidate and that of the 
rest of the tracks and neutral clusters in the event, calcu- 
lated in the T(45) rest frame. The distribution of cos^t 
is sharply peaked near ±1 for combinations drawn from 
jet-like qq pairs and is nearly uniform for the isotropic 
B decays; we require |cos6*t| < 0.9. From Monte Carlo 
simulations of B°B° and B + B~ events, we find evidence 
for a small (1-2%) BB background contribution for the 
channels with rf — > p 7, so we have added a BB compo- 
nent to the fit described below for those channels. 

We use an unbinned, multivariate maximum- likelihood 
fit to extract signal yields and CP-violation parameters. 
We indicate with j the species of event: signal, qq com- 
binatorial background, or BB background. We use four 
discriminating variables: toes, AE, At, and a Fisher dis- 
criminant T [l8j . The Fisher discriminant combines four 
variables: the angles with respect to the beam axis of 
the B momentum and B thrust axis in the T(4S I ) rest 
frame, and the zeroth and second angular moments of 
the energy flow, excluding the B candidate, about the B 
thrust axis 0. 

For each species j and tagging category c, we define a 
total probability density function (PDF) for event i as 

V} tC = Vi (to es 4 ) • Vj (AE*) ■ V, (n ■ V j (Af , oi t ; c) , (2) 

where a\ t is the error on At for event i. With rij defined 
to be the number of events of the species j and fj. c the 



fraction of events of species j for each category c, we write 
the extended likelihood function for all events belonging 
to category c as 

£c = cxp ( - ^2 n j,c) n^s^S'^ig.c ( 3 ) 

j i 
+ n qqfqq,cPqq + n B B I BB ,cPbW ' 

where n J)C is the yield of events of species j found by 
the fitter in category c and N c the number of events of 
category c in the sample. We fix both / s i g , c and f B § c to 
fRty„.r., the values measured with the large -Bfl av sample 
[2fJ. The total likelihood function Cd for decay mode d 
is given as the product over the seven tagging categories. 
Finally, when combining decay modes we form the grand 
likelihood C = Y[ C-d- 

The PDF V s i g (At, o"At;c), for each category c, is the 
convolution of F(At; c) (Eq.QJ with the signal resolution 
function (sum of three Gaussians) determined from the 
-Bflav sample. The other PDF forms are: the sum of two 
Gaussians for "Psig(TOEs) and P s i g (AE); the sum of three 
Gaussians for Vqq (At;c); a conjunction of two Gaus- 
sians with different widths below and above the peak for 
Vj(T) (a small "tail" Gaussian is added for Vqq (J 7 )); a 
linear dependence for Vqq (AE); and for Vqq (to E s) the 
function xV 1 — x 2 exp [—£(1 — x 2 )] , with x = 2toes/v / s- 

For the signal and BB background components we de- 
termine the PDF parameters from simulation. For the 
qq background we use (toes, AE) sideband data to ob- 
tain initial values and ultimately leave them free to vary 
in the final fit. 

We compute the branching fractions and charge asym- 
metry from fits made without At or flavor tagging, ap- 
plied to candidates with ^^^^ and rj' pl combined with 
K + or — > 7r + 7r~. The free parameters in the fit are: 
the signal and qq background yields, the peak position 
and lower and upper width parameters of Vj (J-) for sig- 
nal and qq background, the tail fraction for Vqq (J 7 ), the 
slope of Vqq (AE) and £, the width of the core Gaussian 
of Vgi g (AE), the mean of the core Gaussian of "P S i g (TOEs), 
tibb for B —>■ r/'K, and for charged modes the signal and 
background A c h ■ 

Table |U lists the quantities used to determine the 
branching fraction. Equal production rates of B + B~ and 
B°B° pairs have been assumed. To study biases from the 
likelihood fit, we apply the method to simulated samples 
constructed to contain the signal and background pop- 
ulations expected for data. The resulting yield biases, 
from unmodeled correlations in the signal PDF, are about 
4% for the measurements with 77' — * p°7, and negligible 
for those with 1/^, , . The purity estimate in Table |I] 
is given by the ratio of the signal yield to the effective 
background plus signal, the latter being defined as the 
square of the error on the yield. 

In Fig. ^ w e show projections onto m E s and AE for a 
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TABLE I: Signal yield, purity P(%), reconstruction efficiency 
e(%), daughter branching fraction product, measured branch- 
ing fraction (B) in units of 10 -6 , and Ach- 



Mode 




Yield 


P 




B 


A ch (10" 2 ) 


^77(77)7771 

n ' pl K+ 


K+ 


609 ±28 78 
1347 ± 57 41 
combined 


23 0.175 
26 0.295 


66 ±3 
72 ±3 
69 ±2 


-0.1 ±4.4 
5.5 ±3.6 
3.3 ±2.8 


^77(77)^71 

V 'K° 


K\ _ 

7T — 


198 ± 16 
457 ± 30 
combined 


77 23 0.060 
51 26 0.102 


61 ±5 
73 ±5 
68 ±3 
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FIG. 1: The B candidate jties and AE projections for ri'K + 
(a, b) and if K% (c, d). Points with error bars represent the 
data, the solid line the fit function, and the dashed line its 
background component. 



TABLE II: Results with statistical errors for the B° -> r/'K® 
time-dependent fits. 



Mode 


Signal yield 


S 


C 


(77 ) 7r7r ^ + TT - 


188 ± 15 


0.01 ±0.28 


-0.18 ±0.18 




430 ± 26 


0.44 ±0.19 


-0.30 ±0.13 




54 ±8 


0.79 ± 0.47 


0.11 ±0.35 




44 ±9 


-0.04 ±0.57 


-0.65 ±0.42 




94 ±23 


-0.45 ± 0.68 


0.41 ±0.40 


Combined fit 


804 ± 40 


0.30 ±0.14 


-0.21 ±0.10 




FIG. 2: Projections onto At for 77' K® of the data (points with 
error bars), fit function (solid line), and background function 
(dashed line), for (a) B° and (b) B° tagged events, and (c) 
the asymmetry between B° and B tags. 



subset of the data for which the signal likelihood (com- 
puted without the variable plotted) exceeds a mode- 
dependent threshold that optimizes the sensitivity. 

For the time-dependent analysis, we require \At\ < 
20 ps and cta* < 2.5 ps. We improve the sample size by 
combining the five decay chains listed in Table ITT1 in a 
single fit with 98 free parameters: S, C, signal yields (5), 
r]' pj K° BB background yields (2), continuum background 
yields (5) and fractions (30), background At, m^s, AE, 
T PDF parameters (54). The parameters r and Am^ are 
fixed to world-average values |21| . 

Table [H] gives the yields, S and C, and Fig. the 
At projections and asymmetry of the combined neutral 
modes for events selected as for Fig. ^ 

The major systematic uncertainties affecting the 
branching fraction measurements reflect the imperfect 
knowledge of the rj' branching fractions (3.4%) j2J|, and 
of the reconstruction efficiency (0.8% per charged track, 
1.5% per photon, and 2.1% per K®) estimated from aux- 
iliary studies. We take one-half of the measured yield 
bias (0-2%) as a systematic error. Bias and systematic 
uncertainties for A c h have been estimated from the val- 
ues obtained for the background component in the fit to 
the data. We apply a correction of ±0.016 and assign a 



systematic error of 0.005. 

For the time-dependent measurements, we find approx- 
imately equal (0.01) systematic uncertainties from sev- 
eral sources: variation of the signal PDF shape parame- 
ters within their errors, SVT alignment, position and size 
of the beam spot, BB background, modeling of the sig- 
nal At distribution, and interference between the CKM- 
suppressed b — > ucd amplitude and the favored b — > cud 
amplitude for some tag-side B decays [i^. The Sflav 
sample is used to determine the errors associated with 
the signal At resolutions, tagging efficiencies, and mistag 
rates, and published measurements [2l| for tb and Amj. 
Summing all systematic errors in quadrature, we obtain 
0.02 for S and C. 

In conclusion, we have used samples of about 2000 
rfK + and 800 r/'K® events to measure the branching 
fractions, the time-integrated charge asymmetry and 
the time-dependent asymmetry parameters S and C. 
The measured branching fractions are B(B + — * rj'K + ) = 
(68.9 ± 2.0 ± 3.2) x 10~ 6 and B(B° -> rj'K ) = (67 A ± 
3.3 ± 3.2) x 10~ 6 , and the charge asymmetry is A c h — 
0.033 ± 0.028 ± 0.005. These precise branching frac- 
tion measurements challenge the theoretical understand- 
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ing of these decays [23!]. The measured charge asym- 
metry is consistent with zero, with 90% CL interval 
[—0.012, 0.078], and constrains the amount of possible di- 
rect CP violation in B + — > r)'K + decays. 

The measured time-dependent CP violation param- 
eters in B° -> r/K° are S = 0.30 ± 0.14 ± 0.02 and 
C = -0.21 ± 0.10 ± 0.02. Our result for S differs from 
that measured by BABAR in B° -► J/ipK° by 3.0 
standard deviations; it also represents an improvement 
by a factor 2.4 (1.9) in precision over the published re- 
sults of BABAR 11] (Belle [l2(). All these measurements 
supersede our previous published results 11]. 
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